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Abstract 

The real part of the time-dependent ac susceptibility of the short-range 
Ising spin glass in a transverse field has been investigated at very low tem- 
peratures. We have used the quantum linear response theory and domain 
coarsening ideas of quantum droplet scaling theory. It is found that after a 
temperature quench to a temperature which is less than the spin glass tran- 
sition temperature the ac susceptibility decreases with time elapsed after the 
initial quench towards equilibrium approximately logarithmically. It is shown 
that the transverse field of "tunneling" has nonessential effect to a nonequi- 
librium dynamical properties of the magnetic droplet system. It is found that 
the time dependence of ac susceptibility has a qualitative agreement with 
some experimental results. 
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1 Introduction 



Aging phenomena and nonequilibrium slow dynamics have been investigated during 
last years in many materials with glassy properties such as spin glasses, polimer 
glasses, orientational glasses, simple liquids like glycerol and gels [1-12]. These sys- 
tems are characterized by the existence of a nonequilibrium low temperature phase 
and aging [13-25]. Despite a great progress towards the understanding of nonequi- 
librium dynamics, some problems remain open. One of them is an investigation of a 
very low temperature nonequilibrium dynamics in quantum spin glasses, namely the 
nature of quantum channels of relaxation, the behavior of quantum glassy system 
subjected to periodic driving force, aging at very low temperatures. The natural 
basis for the interpretation of aging is based on coarsening ideas of a slow domain 
growth of a spin-glass type ordered phase [6, 9, 15]. For theoretical studies of 
quantum fluctuations in disordered media there is a variety of techniques including 
replica theory, renormalization group, Monte Carlo simulations, the Schwinger and 
Keldysh closed-time path-integral formalism and others [26-38]. A large attention 
in the last decade was payed to the spin glasses representing a model systems for 
study of nonequilibrium dynamics [39-47]. 

In this paper we investigate the real time nonequilibrium dynamics in d - dimen- 
sional Ising spin glass in a transverse field in terms of droplet model at very low 
temperatures. We calculate the ac susceptibility as a function of the time elapsed 
since a thermal quench. We show that quantum effects alter the nonequilibrium 
dynamics in the spin glass phase at very low temperatures. 

In classical spin glasses in the ac susceptibility measurements the magnetic re- 
sponse of the system to small ac magnetic field after quenching shows aging effects. 
This response depends on its thermal history and the time interval the system has 
been kept at a constant temperature in the glass phase. 

Spin glass materials have nonequilibrium dynamics and provide a measure of 
processes causing the aging. It is assumed that isothermal aging is a coarsening pro- 
cess of domain walls, and the temporal ac susceptibility (real part x' an d imaginary 
part x") a t a given frequency of ac magnetic field oj at time t after the quenching 
scales as [17, 24, 35] 

x"(co,t)-x': q (u J ) 

for |lno;| <C hit if L(t) is proportional logarithm of t; 9 < (d — l)/2. Here L (1/co) 



L(t) 



(1) 



L(t) 



(2) 
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is the typical size of the droplet being polarized by the ac field, L(t) is the typ- 
ical domain size. L may change according to logarithmic growth law (L ~ hit) 
or algebraic law (L ~ t a ) [39, 45]. The relation [L (1/uj) /L(t)] is proportional 
to [In (uj/ujo)/\ii (t/tq)] 1 ^ if the droplet picture is used, ip is a some exponent, 
< ip < d — 1, cuq — Tq 1 , t is a certain microscopic characteristic unit of time. The 
logarithmic growth law (like algebraic growth law) is supported by recent experi- 
ments [17, 45]. The expressions (1) - (2) were found when the relaxation is governed 
by thermal activation over a free-energy barrier B. The barrier for annihilation and 
creation of the droplet excitations are assumed to scale as B ~ AL^; A is a barrier 
energy at T <C T g . The barriers have a broad energy distribution. A droplet with B 
lasts for a time r of order r ~ T exp [B / '(ksT)] where is the Boltzmann constant, 
r is the rate of classical activation over energy barrier B. 

After time t after quenching the typical linear domain size of the system has 
accordingly grown to R(t) ~ [(ksT / A(T))\n (t/ro)] 1 ^. In the ac susceptibility mea- 
surements at angular frequency u, the ac field excites droplets of length scales up 
to L(1/lo) ~ [{ksT I A(T)) |ln(cj/cJo) \\^ ■ Because elapsed time t ^> co^ 1 we have 
L(l/cu) < R(t). These droplets have walls which partly coincide with walls of the 
domain of size R. The presence of such frozen-in domain wall influences the small 
scale droplet excitations. Some droplets which touch it can reduce their excitation 
gap, compared with others in the bulk of domains. In the presence of domain wall 
at a typical distance R from each other D. Fisher and D. Huse have found the free 
energy gap of a droplet of size L in the following form [35] 



t L,R - 7e// 



1 " L<R. (3) 



with effective stiffness constant 7 e // [L/R] =7 1 — c u (L/R) d ~ e , 7 e jj is a function 
of the ratio (L/R). L is a certain microscopic unit of length, a short-distance cutoff; 
c v is a constant independent of time and frequency which happens to be anomalously 
small. So within the droplet picture, aging proceeds by coarsening of domain walls 
as usual phase ordering processes. The domain wall serves as the frozen-in extended 
defect for the droplets and reduces their stiffness constant. It is known that the 
susceptibility is inversely proportional to 7. It was derived, for example, for the real 
part of xi^jt) that [24] 

XV,*) = 4M (l-Cy) (4) 



susceptibility. 



where A7/7 ~ (L(l / uj) / R(t)) [35]. x'eqi 1 ^) is the real part of the equilibrium 
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All aforementioned dynamical processes (droplet excitations, thermal activation 
time and so on) are thermally activated processes. We shall call these processes as 
classical processes in classical regime when (3Tl <C 1. In quantum regime f3Y l ^> 1 
[36, 37], (3 = {ksT)~ l . For quantum droplet model developed in [36-38] L{l/u) ~ 
[(l/a) |ln(r /cj)|] 1/d and L(t) ~ [(l/a) |ln (T t)\] 1/d , where the coefficient a varies 
a little from droplet to droplet. So we suppose the logarithmic growth law of the 
time dependent length scale L(t) of droplet excitations within the frozen-in defect at 
scale R. If we take into account the correction to the true equilibrium which occurs 
due to the finite domain size R we have instead 7 to write 7 e //- Instead of thermal 
relaxation time r for classical process we use in our quantum case the quantum 
tunneling rate for a droplet of linear size L. So the length scale L(t) growth 
is determined not by the thermal activation over the free energy barriers between 
minima but by quantum fluctuations which cause a droplet tunneling through the 
barrier at rates that do not vanish for T — > 0. Thus the quantum fluctuations 
become dominant at very low temperatures. The fraction of droplets which are 
quantum-mechanically active at T — > is proportional to Tl- 

In general, quantum effects could change nonequilibrium dynamics in glassy 
phase at very low temperatures. The understanding aging dynamics in quantum 
system requires a quantum-dynamical approach. 

The paper is organized as follows. In next section we give the definitions and 
main properties of spin glass droplet model. In section 3 we give the linear response 
formalism and general expression for ac magnetic susceptibility. In section 4 we 
present a summary of our results and conclusions. 

2 The model and Hamiltonian 

In this paper we use a phenomenological quantum droplet model of spin glass the- 
ory [35-37] (which does not use the mean-field approximation) in order to describe 
the nonequilibrium behavior of the magnetic dynamical susceptibility at very low 
(but finite) temperatures T (T < IK). We use the quantum Hamitonian for the 
short-range Ising spin glass in a transverse field. This model Hamiltonian may be 
appropriate for experimental systems such as dipolar magnet LiHozYi^Fzi, pro- 
ton glasses, alcali halides with tunneling impurities and other quantum systems [9, 
36-38]. 

We shall be interested in the very low temperatures in the ordered spin glass 
phase and ignore critical effects. The transverse field, applied perpendicular to the 
Ising axis, introduces additional quantum channels of relaxation at that, for example, 
depressing the spin-ordering temperature T g . 
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The droplet model describing the low-dimensional short-range Ising spin glass is 
based on renormalization group arguments [35, 36]. In dimensions above the lower 
critical dimension di (usually in spin glass 2 < d\ < 3) the droplet model finds a 
low temperature spin-glass phase in zero magnetic field. This phase differs essen- 
tially from the spin-glass phase in the mean-field approximation of the Sherrington- 
Kirkpatrick infinite-range spin-glass model [9]. In the droplet model there are only 
two pure thermodynamical states related to each other by a global spin flip. In 
magnetic field there is no phase transition. A droplet is an excited compact cluster 
in an ordered state where all the spins are inverted. The natural scaling ansatz for 
droplet free energy (which are considered to be independent random variables) 
is €l ~ L e ,L > ((T); ( is the correlation length, L is the length scale of droplet 
and 9 is the zero temperature thermal exponent. One droplet consists of order L d 
spins. Below di, 9 < 0; above di one has 9 > 0. The droplet excitations have a 
broad distribution of their free energies with the probability distribution Ph{^L)deh 
of droplet free energies at scale L for large L in a scaling form [4, 6] 

P de L )de L = [-^) ,L^oc. (5) 

It is assumed that P L (x -> 0) > 0, V L (0) - V L {x) ~ x* at x -> 0; < 1. 7(T) is the 

stiffness constant of domain wall on the boundary of the droplet which is of order 

i 

of characteristic exchange X = (xf^j 2 at T = and vanishes for T > T g . 

The Hamiltonian of the d - dimensional quantum Ising spin glass in a transverse 
field is given by 

H = -Y t ZaS?S' i -TY t S! (6) 

i,j i 

where Si are the Pauli matrices for a spin at site %. T is the strength of the transverse 

field. The sum in (6) is performed over nearest neighbors. The interactions 

l. 

are independent random variables of mean zero and variance X = (xf^j 2 . This 
model shows, for example, the physics of proton glasses, mixed betaine phosphate- 
phosphite [9]. The transverse field may be interpreted as the frequency of proton 
tunneling. The quantum spin glass transition in a dilute dipole coupled magnet 
LiHo :r Y 1 _ :E F4 [9] also is described by this model Hamiltonian. The properties of 
model (6) in mean-field approximation have been studied in many papers [6, 9]. 
It was found that there is a high-temperature critical behavior at the temperature 
T c (r) ~ X, T X and low-temperature critical behavior with the zero-temperature 
critical point T C (X) = 0, r c (0) = X, where r c is the critical value of T below which 
the spin-glass phase can exist. We suppose as in our early papers [37, 38] that 
quantum system with the Hamiltonian (6) has a true glass phase transition. 
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We shall use the quantum model Hamiltonian (6) of the d - dimensional Ising 
spin glass in the transverse field. If T = the Hamiltonian (6) is the d - dimen- 
sional classical Ising spin glass. One can use the Suzuki - Trotter formalism [36] 
to show that the d - dimensional quantum mechanical system (6) is equivalent to a 
classical statistical mechanical system in (d + 1) - dimensions (the extra dimension 
is imaginary time) with the classical Hamiltonian 

H d = At X X ZijSi,jSj,k ~ E X S itk S itk+1 (7) 

k=l (ij) fc=l i 

where the S^k = ±1 are the classical Ising spins, representing the ^-component of the 
quantum spins at site i and imaginary time r = A; At, the indices i and j running over 
sites of the d - dimensional lattice. The imaginary time direction has been divided 
into L T time slices of width At. The calculation gives that exp (1 F ) = tanh (AiT), 
AtL t = (ksT)" 1 , here and throughout the paper we use units where H — 1. Using 
properties of the Hamiltonian (7) of the (d + 1) - dimensional classical model at 
temperature T = M. J. Thill and D. A. Huse [36] have assumed that in each of the 
two ordered states (for r < T c and T < T g ) for enough low T and appropriate range 
of length scales L the quantum Hamiltonian (6) can be represented as independent 
quantum two-level systems (low energy droplets) with the Hamiltonian 

w = ^EE( c ^ + r ^J ( 8 ) 

L Dl 

where S z Dl and Sf, are the Pauli matrices representing the two states of the droplet; 
the sum is over all droplets D L at length scale L and over all length scales L, and 

roc dL 
l JL ° 

where Lq is a short-distance cutoff. eo L is the droplet energy which is independent 
random variable. The droplet length scale L is more or of order of the correlation 
length (. 

r L = r ex P [- ( 7L d ] (io) 

is the tunneling rate for a droplet of linear size L, T is the microscopic tunneling 
rate; a is the surface tension for the interface between the two droplet states, a 
is approximately the same for all droplet. We will assume Y ^ is the same for all 
droplets of scale L [36] . The Hamiltonian of a single droplet is the 2x2 matrix 



~ dL 



( r \ 



y T L -6 Dl j 



(11) 
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with eigenvalues E± = ±y e| + Y\. E = 2\E±\ is the energy difference between the 
two eigenvalues. Note the Hamiltonian (9) is similar to the Hamiltonian of two-level 
system in real glass [32]. 

In the quantum droplet model of M. J. Thill and D. A. Huse [1] the relative 
reduction of the Edwards-Anderson order parameter q B A(T) from its zero-T value 
q E A(0) for 6 > and L*(T) > ( is given by 

q EA (T) k B T k B T a e / d 

q EA {T = 0) ~ 7 L*(T) 7 ^ d (F /(k B T)) 1 ' 1 ' 

where the crossover length scale is 

L*(T)= (-\n-%Y . (13) 



.a k B T 

For droplets with L L*(T) and T L » k B T the energy \Je 2 L + Y 2 L is always 
greater than k B T and thermal fluctuations are therefore not essential at temperature 
T. These droplets behave quantum-mechanically. The larger droplets (L 3> L*(T)) 
have T L <ti k B T and behave classically. Those large droplets (e_L < k B T, T L < k B T) 
are thermally active and reduce q B A- For 6 > this reduction is dominated by 
smallest droplets. Below d t (9 < 0) the thermally-excited droplets disorder the 
system reducing q B A to zero at any temperature. 

There is complicated dynamical-to-quantum crossover depending on tempera- 
ture, frequency of ac field and length scale L. According to [36] the crossover 
dynamical length is determined from condition T^ 1 = t, i. e. 

L' im (T)~(^k B T)^. (14) 



The system behaves presumably classically or quantum-mechanically when the dom- 
inant length scale L is above or below L* dyn for fixed frequency u. When the droplets 
behave quantum-mechanically they have a characteristic rate - the Rabi frequency 
(is of order T B )- 

In this paper we neglect droplet-droplet interaction and droplet-lattice one. 

3 Linear response 

We consider the time-dependent Hamiltonian 7i of the quantum system in the form 
[48] 

H = H + H'(t) =H - Ah{t). (15) 
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where Hq is nonperturbated Hamiltonian of system, H ^ Ti-(t) and suppose that 
the external perturbation 7i'(t) is in some sense small. 7io describes the equilibrium 
system. Here the A is the linear operator through which the external time varying 
force h(t) couples to the system. 

We use the quantum-mechanical calculation for the system dynamical response 
AB(t) = (B(t)) — (B)o to the force h(t) in terms of the time-evolution operator 
U(t,t'); B(t) is the Heisenberg operator, B(t) = U^(t,t')B(t )U(t,t'), (B) is the 
equilibrium expectation value of B. It is necessary to make approximation for U (t, t') 
carrying out the well-known perturbation expansion of U (t, t') through first order 
in applied force in the following form 

u(t,t') ~ u (t,t') ji - ^ jfV# t (*i,Ow , (fi)#o(*i,o} (16) 

where U (t,t') = exp — | (t — t')H , the sign f means the conjugate value. 

We take into account the first- and second-order linear response functions for 
the linear response functional of the form [48] 

$(t, t'- 1 ) = ±.([A(t ), B(t,t')}) (17) 

where (. . .)o means the thermal average with the density matrix p = p(to), where to 
is the time moment when the perturbating field is turned on, p = (Tr exp[— /37i]) _1 
x exp[-f3H]. 

Now we apply the aforementioned dynamical response relations to our magnetic 
droplet system. Then the response (B(t)) represents the induced magnetization of 
droplet system M(t). (B) is the equilibrium magnetization M . If in z-direction a 
small magnetic oscillating field h(t) = hexp[iut] is applied, the induced magnetiza- 
tion (in ^-direction) of the sample is 

M(t) = {x'(lu) cos(cut) + x"H sm(cut)}hV, (18) 

where x{°°) = x'{ u ) ~ i^'i 00 ) an d V is the sample volume, Xzz = X\ h and u is the 
amplitude and the angular frequency of ac field. 

In order to observe a history dependence and aging in spin glass the sample is 
quenched infinitely fast at zero dc magnetic field from temperature T ^> T g to the 
temperature 7\ < T g which is reached at time t = 0. At this moment a very small 
external magnetic oscillating field h(t) is applied to measure the ac susceptibility of 
the sample. The evolution continues in isothermal conditions, Xac is measured as a 
function of the time t elapsed since the sample reached the temperature 7\ at fixed 
frequency uj. 
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If the field changes are small the magnetization M(t) responses linearly to a 
magnetic field h(t) applied at once when the quench was finished, i. e. at time 
moment t — 0. The system is probed at the time t after quench end (the "age"). 
Using linear response theory the magnetization of droplet system is [25] 

M(t)-M = f x{t,t 1 )h{t 1 )dt 1 = f x(t,t-t')h(t-t')dt' (19) 
Jo Jo 

where x(^ t — t') is the magnetic dynamical susceptibility (the dynamical response 
functional including the first- and second-order linear response functions). xit^t — t') 
determines the magnetic response at time t to a unit magnetic field impulse at 
time (t — t'). The nonequilibrium processes are probed by the low-frequency ac 
susceptibility measurements. The frequency dependent ac susceptibility is measured 
by applied ac magnetic field h(t) at time t — 0. Then x(w,t) may be found by the 
Fourier transform of the magnetization over time segment t m (t m ~ 2ti/uj) centered 
around t [25, 39] 

1 rt+— rt" 

X (u, = 7-// dt"e-^" / dt' X (t", t" - t')e^"-^ (20) 

tm Jt-if- JO 

If the magnetic response function varies little over time segment t m the susceptibility 
x{w,t) will be equal to [25] 

X (u,t)= f 'dt'x(t,t-t')c-^' (21) 
Jo 

The in-phase component of the ac susceptibility x'(w,t) can be written 



X '(u,t)=Re\ [* dt' X (t,t-t'] 
Jo 



e -iut 



The out-of-phase component of the ac susceptibility is given by 



X "(u;,t) = lm\f t dt'x(t,t-t') 
Jo 



(22) 



(23) 



We consider the behavior of the magnetic droplet system described by the Hamil- 
tonian H under ac field h(t) in quantum regime (T L ^> k B T) when the droplet 
excitation energy \Je 2 L + Y 2 L is always greater than k B T, i. e. the temperature is too 
small to thermally excite the droplets which behave classically and thermally active 
(classical regime). In our calculations below we assume that 9 > (d > d{). 

Following to aforementioned quantum droplet theory and domain growth ideas 
we have the following expression for susceptibility xd l (xd l is proportional to the 
contribution of a single droplet to the ac susceptibility) up to some factor ~ g| A L 2d 

XD L (w,t) - Xd l {u = 0) ~ -tanh (J3a L /2) sin 2 ip ^ 2 

a L — lo 
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— /itanh (/3ol/2) cos y2 sin 2 y9 



(7a 2 L - I0u 2 ) cos{out) 



3a,L sm(ut) sin(a^) + 6ou cos(ut) cos(a£,i) cos(ar,£) 
Ld(a 2 L — Auu 2 ) a 2 L — uo 2 



+i 



3(a| + 2cj 2 ) sin(cut) 
*| -^ 2 )(a| -4a> 2 ) 



3ol cos(a;t) sin(ait) — 6cj sin(cjt) cos(aLt) 3 sin(ait) , , 

+ ~ 772 ~77T f ' 




u;(a2 — 4c; 2 ) a;(a| — w 2 



where = ye| + r|, siny? = T L /a L , cosip = t^ja L , Xd l (u = 0) is the static 
susceptibility of droplet D^. The expression (24) was obtained with condition out > 1 
(and Tl < uj) because this condition is used for observing nonstationary dynamics 
in Xac measurements [12]. If out > 1, i. e. for low frequencies, the aging part 
usually strongly contributes. For high frequencies, out — > oo, the aging term does 
not contribute. Now we have to average the susceptibility (24) over droplet energies 
ei and over droplet length scales L. In order to average over droplet energies we use 
the distribution of free energies Pl{cl) (5). In this distribution we assume = 0. We 



approximate in (24) tanh(ye| + T L /2) ~ 1 and integrate like in papers [36, 37, 38]. 
As the result of the droplet energy averaging we received the averaged contribution 
of all droplets of the system to the real part of susceptibility in the following form 



X L (u,t)-x L {uJ = 0) In \ -In- 



^uyfd* - r 2 \ u- ^/uj 2 - r 2 ' uj Pl - sju 2 - vlh 



Y 2 h r /3 3 

^jjj, [sin(c^) ^-si (b L t) + si(t(6 L - w)) + si(*(6 L + ou)) + -si (f (& L - 2oo)) 



3 \ / 3 3 

+ -si (t(b L + 2cj))J + cos(wt) I --ci (t(& £ - 2uo)) + -ci (t(& £ + 2uo)) - In 



&l + ou 



ou 



3 &£, + 2oj 5uj 3uj . 
--In ci (t(b L - uj)) + ci (£(& £ + a;)) + — cos(b L t) - 3utsi (b L t) 

4 b L -2uj b L b L 



cos(b L t) — 2utsi (b L t) 

bL 



(25) 

where XzX^ = 0) * s the static susceptibility of the system of droplets of size L, 



bj_, = \Jp 2 L + r 2 , p£ = 2ou + r L , si(a) is the sine integral and ci(a) is the cosine 
integral. Further we average expression (26) over length scales L. While integrating 
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over L we see that the real part of the susceptibility is dominated by droplets of 
length scale L\ ~ [(l/a)\n(To/uj)] 1 ^ d . L\ is the natural length scale of the problem 
when ~ uj and the lower limit of integration over L. The upper limit of integration 
we took as L2 ~ [(1/<t) ln^ro)] 1 ^. Hence we obtain after averaging over L the 
following expression for the real part of the ac susceptibility of droplet system 



xW)-xV = o) 



r 2 
1 



(ln^-2)(2a)^- 1 G 
1 



e 



m* 

UJ 



+ ^(21n^-l)(4a)^- 1 G 







In* 

UJ 



9 

T 6 



m* 



UJ 



9 

r 2 



UJ 



UJ 



9 

r 4 



>n* 

UJ 



r^/i f f 1 3 17 3 

H — ^-r ( ( — (— sin(3u;t) — — sin(2ujt) — - sin(c<jt) cos(2cc;t)) + cos(o>t) 
'yuj Wujt 16 12 4 



x(5-ln3-|ln^)) (-d- 1 (2,)l(G 



•^,2|ln(tr ) 



G 



(Tcos(cjt) ^-rf" 1 (2a)^ 1 (G 1 - ^ , 2 |ln(*r 



G 



,2 



In 



r 



UJ 



(26) 



where G[a, x] is the incomplete gamma function. 

If we use the asymptotic representation for the incomplete gamma function for 
large values of the its second argument we obtain for the difference AG of two 
incomplete gamma functions 



AG 



1 - 



ln^- 1 ^) 



1+ 



UJ 

f 2 ' 
1 



(27) 



(ujty v iu(*r ) ; 

We observe a some similarity with expression (5.7) (in brackets) in [35] in con- 
formity with our quantum regime. 

At derivation of (26) we made the following approximations: 1) si(6z,t) ~ —{Pht)^ 1 
x cosGM); 2) ci(M) ^ (prf)- 1 sm(p L t); 3) b L ~ 2uj + T L . 

We also made the numerical calculation of the average over L of the expres- 
sion (25) without these approximations and obtained almost the same curves for 
susceptibility as ones obtained from (26) (the parameter values we give below). 
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The susceptibility x'(u;,t) depends on many parameters of the droplet system 
and the external ac magnetic field, for example, it depends on the kind of distri- 
bution function Pl{cl), on the droplet microscopic tunneling rate Tq, on time t, on 
frequency and amplitude of ac magnetic field, on the droplet features. 

The expression (26) consists of the time-independent terms, the terms which 
describe the simple oscillations with frequency 00 and the terms which depend on 
time t and determine nonstationary nonequilibrium dynamics of droplet system. So 
the real part of the ac susceptibility can be represented approximately as a sum of 
the stationary part (x'st) an d nonstationary part (x'nst)'- 



Let us take for numerical calculation of expression (26) the following values of 
parameters: d = 3, 9 = 0.5, 7 = KT 15 , T = 10 8 , 10 10 , 10 12 , h = 1, a = KT 15 , 
* = 0-5-100, co = 0.05, 0.1. 

In Fig. 1 we show the t - dependence of the real part of ac susceptibility of 
the droplet system x'nst- Here out is comparable or more than unity, so one may 
observe nonstationary dynamics and the aging regime [12]. In Fig. la it is observed 
the slow dynamics at r = 10 s , to = 0.05, 0.1 and t = 4- 25. At first the curve 
grows (very quickly) up to some finite value and then falls down. At more long 
times the ac susceptibility shows a decay in queue. In Fig. lb and Fig. lc the t - 
dependence of x'nst( u ^ t) * s shown for more long times: b) t — 04-50; c) t = 04-100. 
The time interval covers two decades of the elapsed time t. We see that with 
frequency increasing the susceptibility magnitude decreases and the slow dynamics 
is suppressed at more high frequency. The ac susceptibility tends more rapidly to 
zero at more long times. In Figures 1 we see the essential influence of frequency 00 
and elapsed time t on the susceptibility. 

In Fig. 2 we give the time dependence of x'nst( u ^) a ^ different values of fixed 
quantum parameter r (the microscopic tunneling rate r = 10 8 , 10 10 , 10 12 ). We 
observe the small effect of quantum parameter r on susceptibility x'nst(, u i With 
increasing of r the magnitude of the susceptibility slightly decreases at small times, 
i. e. the quantum fluctuations in some sense decrease the susceptibility in quantum 
regime in spin glass phase. 

We may compare our data with experimental ones [18, 24, 25] only very ap- 
proximately because in the x'i^i t) ~ experiments and in our paper the qualitatively 
different spin-glass systems are considered. We observe qualitatively similar dynam- 
ical behavior of %'(cj,t) in the range of the small times elapsed since the quench. 
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4 Discussion and conclusion 



In this paper we have investigated the low temperature nonequilibrium dynamical 
behavior of the magnetic ac susceptibility in d - dimensional Ising spin glass with 
short-range interactions between spins in a transverse field in terms of the phe- 
nomenological droplet model. The real part of the low-frequency ac susceptibility 
x'(ou,t) as a function of time t (which the system waits from quench time moment 
till the time moment of measurement) and the frequency ui of external ac magnetic 
field is calculated. We display the nonequilibrium dynamics for different values of 
low frequency of ac field at constant temperature in spin glass phase. We can see 
the oscillations of x' t) from expression (26). The real part of ac magnetic suscep- 
tibility x'(w,t) of droplet system at very low temperatures (quantum regime) has 
two-time regions in which the time evolution is of a different nature. At short times 
t we observe nonequilibrium dynamical behavior - the time decay of x'(w,t) at low 
frequency and at constant temperature T (Tl ^ ksT in quantum regime). If the 
frequency of ac field increases the nonequilibrium dynamics is suppresed. So, the 
response of droplet system to external perturbating field depends on the thermal 
history. 

In [44] it was shown that the behavior of the response function R(t, t w ) demon- 
strates the existence of the stationary and aging regimes in quantum systems. The 
theoretical curve (Fig. 2 in [44]) R(t,t w ) was given as function of r (r = t — t w ) 
if r G [0, 50] and t w = 2.5, 5, 10, 20 and 40. For r < r char [r char is a some char- 
acteristic time) it was found the stationary regime, if r > r char the dynamics is 
nonstationary. In [18, 24, 25] the x'(^,t) - experimental data are represented. In 
aging regime the slow monotonous decay of x'(w,t) was observed. In our quantum 
system at very low T we can not find good agreement with these data. 

Besides we have shown that the quantum fluctuations have slight influence to 
the dynamical susceptibility of droplet system at very low temperatures. Their 
increase leads to the small decrease of the susceptibility magnitude. Like in [43] 
it is shown that in the aging regime of quantum spin glasses of rotors all terms 
in the dynamical equations governing the time evolution of the spin response and 
correlation function that arise solely from quantum effects are irrelevant at long 
times. The quantum effects enter only through the renormalization of the dynamical 
equations parameters [43]. In paper [44] it is shown that quantum fluctuations in 
quantum glassy systems depress the phase transition temperature, in a glassy phase 
aging survives the quantum fluctuations and the quantum fluctuation - dissipation 
theorem is modified due to quantum fluctuations. 

Thus, we have found the nonequilibrium low-frequency dynamical behavior of 
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the ac susceptibility at constant temperature for droplet system in quantum regime 
(r^ ^> fcfiT) in spin glass phase. For small times t the magnetic ac susceptibility de- 
pends on both frequency of ac magnetic field u and on time t elapsed since the sample 
reached the temperature T\ (Ti < T g ). We do not find slow continues decrease of 
the amplitude of x'{^,t) as a function of the time t at long times. Aging is more 
visible in the out-of-phase component x"(u;,t) of the magnetic susceptibility which 
determines a dissipation [12]. A qualitative agreement with x'(u;, t)-measurements 
[18, 24, 25] in spin-glass systems was found for small times and low frequencies. 
This work is partially supported by the RBRF under Grant 01-02-16368. 
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Figure captions 



Fig. 1. The in-phase susceptibility x'(cj,t) (nonstationary part) as a function of 
time t for quantum parameter r = 10 8 and fixed values of frequency uj = 0.05, 0.1: 

a) uj = 0.05, 0.1; t = -=- 25; 

b) uj = 0.05, 0.1; t = 0^50; 

c) cu = 0.05, 0.1; t = 0^100; 

Fig. 2. The in-phase susceptibility x'(c<j,t) (nonstationary part) as a function 
of time t for frequency ou — 0.1 and for three values of quantum parameter Tq = 
10 8 , T = 10 10 , T = 10 12 . 
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